Indications of a spatial variation of the fine structure constant 
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We previously reported Keck telescope observations suggesting a smaller value of the fine structure 
constant, a, at high redshift. New Very Large Telescope (VLT) data, probing a different direction 
in the universe, shows an inverse evolution; a increases at high redshift. Although the pattern could 
be due to as yet undetected systematic effects, with the systematics as presently understood the 
combined dataset fits a spatial dipole, significant at the 4.2a level, in the direction right ascension 
17.5±0.9 hours, declination — 58±9 degrees. The independent VLT and Keck samples give consistent 
dipole directions and amplitudes, as do high and low redshift samples. A search for systematics, 
using observations duplicated at both telescopes, reveals none so far which emulate this result. 

PACS numbers: 06.20.Jr, 95.30.Dr, 95.30.Sf, 98.62.Ra, 98.80.-k, 98.80.Es, 98.80.Jk 



Quasar spectroscopy as a test of fundamental 
physics. — The vast light-travel times to distant quasars 
allow us to probe physics at high redshift. The relative 
wavenumbers, uj z , of atomic transitions detected at red- 
shift z — Xobs/^iab — 1, can be compared with laboratory 
values, loq, via the relationship uj z = u>o+Q (a 2 — a 2 ,) /a 2 , 
where the coefficient Q measures the sensitivity of a given 
transition to a change in a. The variation in both magni- 
tude and sign of Q for different transitions is a significant 
advantage of the Many Multiplet method [U [2] , helping 
to combat potential systematics. 

The first application of this method, 30 measurements 
of Aa/a = (a z — ao) /uq, indicated a smaller a at high 
redshift at the 3er significance level. By 2004 we had made 
143 measurements of a covering a wide redshift range, 
using further data from the Keck telescope obtained by 
3 separate groups, supporting our earlier findings, that 
towards that general direction in the universe at least, a 
may have been smaller at high redshift, at the 5<r level [3]- 
[5]. The constant factor at that point was (undesirably) 
the telescope and spectrograph. 

New data from the VLT. — We have now analysed 
a large dataset from a different observatory, the VLT. 
Full details and searches for systematic errors will be 
given elsewhere [6j [7]. Here we summarize the evidence 
for spatial variation in a emerging from the combined 
Keck+VLT samples. Quasar spectra, obtained from the 
ESO Science Archive, were selected, prioritising primar- 
ily by expected signal to noise but with some preference 
given to higher redshift objects and to objects giving 
more extensive sky coverage. The ESO midas pipeline 
was used for the first data reduction step, including wave- 
length calibration, although enhancements were made to 
derive a more robust and accurate wavelength solution 
from an improved selection of thorium-argon calibration 
lamp emission lines [8] . Echelle spectral orders from sev- 
eral exposures of a given quasar were combined using 
tjves_popler [5]. A total of 60 quasar spectra from the 



VLT have been used for the present work, yielding 153 
absorption systems. Absorption systems were identified 
via a careful visual search of each spectrum, using rd- 
GEN [TU] , scanning for commonly detected transitions at 
the same redshift, hence aligned in velocity coordinates. 
Several transition matches were required for acceptance 
and, given the high spectral resolution, chance matches 
were eliminated. 

Absorption system modelling. — As in our previous 
studies, vpfit was used to model the profiles in each ab- 
sorption system with some enhancements, described 
in [BJ. A comprehensive list of the transitions used, their 
laboratory wavelengths, oscillator strengths, and Q coef- 
ficients are compiled in [BJ . 

The following general procedures were adhered to: (i) 
For each absorption system, physically related param- 
eters (redshifts and ^-parameters) are tied, in order to 
minimise the required number of free parameters and 
derive the strongest possible constraints on line posi- 
tions, and hence Aa/a. (ii) Parameters were tied only 
for species with similar ionisation potentials, to min- 
imise possible introduction of random effects on a, mim- 
icked by spatial (and hence velocity) segregation effects; 
(Hi) Line broadening is typically dominated by turbulent 
rather than thermal motion. Both limiting-case models 
were applied and Aa/a determined for each. The final 
Aa/a was derived from a likelihood- weighted average; 

(iv) Where appropriate and available, isotopic shifts and 
hypcrfine structure are included in the fitting procedure; 

(v) Velocity structures were determined initially choos- 
ing the strongest unsaturated transitions in each system. 
Normalised residuals across each transition fitted were 
examined and the fit progressively refined with the intro- 
duction of each additional transition to the fit; (vi) Tran- 
sitions falling in spectral regions contaminated by telluric 
features or atmospheric absorption were discarded. Any 
data regions contaminated by cosmic rays, faulty CCD 
pixels, or any other unidentified noise effects, were also 
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discarded; (vii) A few gravitational lenses were identified 
by being difficult or impossible to model successfuly. The 
non-point source quasar image and the resultant complex 
line-of-sight geometry can significantly alter apparent rel- 
ative line strengths. These systems were discarded; (viii) 
In all cases we derived the final model without solving for 
Aa/a. The introduction of Aa/a as an additional free 
parameter was only done once the profile velocity struc- 
ture had been finalised, eliminating any possible bias to- 
wards a 'preferred' Aa/a. One potential consequence 
of this approach might conceivably be a small bias on 
Aa/a towards zero, should some 'fitting-away' of Aa/a 
occur by column density adjustments or velocity struc- 
ture decisions. The reverse is not true, i.e. it cannot bias 
towards a non-zero Aa/a. For details of all the points 
above see [B] . 

vpfit [IT] minimises x 2 simultaneously over all species. 
Whilst the strongest components may appear in all 
species, weaker components can sometimes fall below the 
detection threshold and hence are excluded, such that 
a component which appears in Mgll, for example, does 
not appear in Fell. There is no solution to this (known) 
problem but its effect merely adds an additional random 
scatter on Aa/a for an ensemble of observations. 

Spatially dependent a. — An initial inspection of Aa/a 
vs redshift for the new VLT dataset reveals a redshift 
trend, opposite in sign compared to the earlier Keck 
data. Splitting each sample at z = 1.8, our 2004 Keck 
sample [5] gave (Aa/a),<i.g = -0.54 ±0.12 x 10 -5 and 
(Aa/a) z>1 . s = -0.74 ± 0.17 x 10~ 5 . The present VLT 
sample, discussed in detail in [6], gives (Aa/ a) z< \.% = 
-0.06 ± 0.16 x 10~ 5 and (Aa/aj z>1 . 8 = +0.61 ± 0.20 x 
10~ 5 . Errors here and throughout this paper are la es- 
timates. Our VLT result above for z < 1.8 agrees with 
the VLT data presented in |12) . 

Errors on individual Aa/a values for our VLT sam- 
ple are a 2 tot = a 2 tat + a 2 rand , where a 2 rand was derived 
empirically by fitting a constant Aa/a to the sample, 
i.e. monopole-only, using a modification of the Least 
Trimmed Squares (LTS) method, where only 85% of 
data, those points with the smallest squared residuals, 
are fitted. o- ran d was assumed constant for all absorbers 
and found to be s» 0.9 x 10 -5 , showing that the scat- 
ter in the VLT Aa/a is greater than expected on the 
basis of statistical-errors alone. Errors on Aa/a for the 
Keck sample are discussed in [4] , although we derive a 
new estimate of o ra nd = 1.74 for the Keck points using 
the LTS method, again relative to a monopole-only fit 
to the Keck sample. 

The Keck (Mauna Kea, Hawaii) and VLT (Paranal, 
Chile) locations on Earth are separated by 45° in latitude 
and hence, on average, observe different directions on the 
sky. The (Aa/a) results above suggest exploring a simple 
spatial dependence using the combined dataset. 

The Keck sample we use is as presented in [5] with 
minor modifications: 3 points were removed. 2 had been 



included erroneously (from a spectrum known to have 
calibration problems) and 1 further point was clipped, 
having a residual greater than 3er against a modified LTS 
fit to the Keck data. 

We fit 3 different models to the 3 datasets (i.e. 
Keck, VLT and combined samples). Initially we try a 
dipole -I- monopole model, Aa/a— AcosO + m, where 
m allows an offset from the terrestrial value, O is the 
angle on the sky between quasar sightline and best-fit 
dipole position, and A is the dipole amplitude. Noting 
the theoretical interpretation of the monopole term is 
unclear, we fit a second model, without the monpole, 
Aa/a= A cos©. Thirdly, in order to explore a possible 
spatial gradient in a, we assign a distance to each Aa/a 
measurement of r(z) = ct(z) where c is the speed of light 
and t(z) is the look-back time at redshift z. The model 
is then Aa/a= ArcosQ. 

To estimate the dipole significance we bootstrap the 
sample, repeatedly randomising the association between 
Aa/a and the absorption system location in space (i.e. 
quasar sightline and absorption redshift). A dipole is fit- 
ted and its \ 2 derived at each realisation, to obtain a 
X probability distribution. This gives the probability 
of fitting a dipole to the data and obtaining a value of 
X 2 less than or equal to that observed for the real sam- 
ple by chance alone, i.e. the statistical significance of 
a dipole model compared to a monopole, and hence an 
uncertainty estimate for the dipole amplitude, A. 

All 3 models give a detection significance in the range 
4.1 — 4.2cr and the best-fit parameters and associated 
errors (given in the figure captions) vary only slightly. 
Figure 1 illustrates an all-sky map for the best-fit no- 
monopole model, using equatorial co-ordinates. Approx- 
imate la error contours are derived from the covariance 
matrix. Figure 2 illustrates the Aa/a binned data and 
the best-fit dipole+monopole model. Figure 3 illustrates 
Aa/a vs look-back time distance projected onto the 
dipole axis, rcosO, using the best-fit dipole parameters 
for this model. This model seems to represent the data 
reasonably well and Aa/a appears distance-dependent, 
the correlation being significant at the 4.2<r level. 

An alternative to the LTS method described above, to 
allow for any unknown additional contribution to the er- 
rors on individual Aa/a measurements, one can assume 
a\ ot = a 2 tat and iteratively trim the sample during model 
fitting. This provides a further test of whether the ap- 
parent spatial gradient in a is dominated by a subset of 
the data, perhaps more prone to some unknown system- 
atic than the remainder. Adopting a 2 ot = a 2 tat will tend 
to result in higher significance levels. Figure 4 illustrates 
this test and shows that the apparent dipole seems robust 
to data trimming. 

Empirical test for systematics. — One potential sys- 
tematic in the data could arise if there were slight me- 
chanical mis-alignments of the slits for the 2 arms of the 
UVES spectrograph on the VLT. This could cause wave- 
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FIG. 1. All-sky plot in equatorial coordinates showing the in- 
dependent Keck (green, leftmost) and VLT (blue, rightmost) 
best-fit dipoles, and the combined sample (red, centre), for the 
dipole model, Aa/a= AcosQ, with A = (1.02L0.21) x ICT 5 . 
Approximate la confidence contours are from the covariance 
matrix. The best-fit dipole is at right ascension 17.4 ± 0.9 
hours, declination —58 ± 9 degrees and is statistically pre- 
ferred over a monopole-only model at the 4.1a level. For this 
model, a bootstrap analysis shows the chance-probability of 
the dipole aligments being as good or closer than observed 
is 6%. For a dipole+monopole model this increases to 14%. 
The cosmic microwave background dipole and antipole are 
illustrated for comparison. 
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FIG. 2. Act/ a for the combined Keck and VLT data vs angle 
from the best-fit dipole position (best-fit parameters given 
in Figure 1 caption). Dashed lines illustrate ±1<t errors. For 
a discussion on the monopole term, see [6j . 



length shifts between spectral features falling in the blue 
and red arms. However, this specific effect appears to be 
substantially smaller than required to explain values of 
Aa/a~ 10 -5 seen in the present work [T3"] . 

A more subtle but related effect may be slight off- 
centre placement of the quasar image in the spectrograph 
slit, by different amounts for different exposures, at dif- 
ferent wavelength settings. This may apply to either or 
both Keck and VLT spectra. Since spectrograph slit il- 
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FIG. 3. Act/ a vs Arcos0 showing an apparent gradient in 
a along the best-fit dipole. The best-fit direction is at right 
ascension 17.5 ± 0.9 hours, declination —58 ± 9 degrees, for 
which A = (1.1 ± 0.25) x 10~ 6 GLyr -1 . A spatial gradient 
is statistically preferred over a monopole-only model at the 
4.2<r level. A cosmology with parameters (Ho,flM,^A) = 
(70.5,0.2736,0.726) was usedQH]. 
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FIG. 4. As an alternative to increasing Act/a error bars, 
to account for the additional scatter in the data as described 
in the text, we instead use a\ ot = a^tat and iteratively clip 
the most deviant Act/a value, fitting Aa/a— Arcos0. Ap- 
proximately 60% of the data must be discarded before the 
significance drops below 3a showing the dipole signal is not 
due to a small subset of the data. The solid (pink) line at 
the bottom of the graph shows the dipole amplitude in units 
of 10~ 6 Glyr" 1 . The dotted (blue) line at the bottom of the 
graph shows xi an( A the vertical dashed (red) line illustrates 
xt = 1 when ~ 8% of the data has been trimmed, at which 
point the significance is ~ 5.5<r 



luminations are different for quasar (point source) and 
ThAr calibration lamp (uniform illumination), the sub- 
sequent combination of individual exposures to form a 
1-dimensional spectrum may then contain relative veloc- 
ity shifts between spectral segments coming from differ- 
ent exposures. This effect will exist in our data at some 
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level and it is clearly important to know the impact on 
an ensemble of measurements of a. 

Fortunately, 6 quasars in our sample have both Keck 
and VLT spectra, allowing a direct and empirical check 
on the effect above, and indeed any other systematics 
which produce relative velocity shifts along the spectrum. 
To do this we selected small spectral segments, each a few 
A wide, flanked by unabsorbed continuum flux from the 
quasar, and fitted Voigt profiles using vpfit, but adding 
an additional free parameter allowing a velocity shift be- 
tween the Keck and VLT segments, Sv(X b s )i, where X i, s 
is the observed wavelength and i refers to the i th quasar. 
All available absorption lines in the 6 spectra were used, 
including both Lyman-a forest lines and heavy element 
lines but excluding telluric features. In this way we can 
map any effective relative distortions in the calibrations 
between each pair of spectra. A total of 694 measure- 
ments were used from the 6 pairs of spectra over the 
observed wavelength range 3506 < A < 8945A. 

We formed a composite function Sv(X b s ) after first 
normalising (Sv(X b s )i) = for each i to remove any po- 
tential small constant velocity offsets from each spectrum 
(expected from off-centering of the quasar in the spectro- 
graph slit), which cannot influence a. 

Finally, we fit the composite 6v(X b s ) with a linear 
function f(Sv) = aX i. s + b where a — (—7 ± 14) x 
10 _5 kms- 1 A- 1 , b = 0.38±0.71 kms" 1 . The final f{Sv) 
thus shows a weak (but statistically insignificant) veloc- 
ity drift, and provides an empirical transformation be- 
tween the Keck and VLT wavelength scales. For each 
VLT quasar absorption system, we modify the input 
laboratory wavelengths used in the Voigt profile fitting 
procedure Xi ab to X' lab = X lab + AA /a6 where AA/ a& = 
Xiab 5v (X bs) / c, and finally use the X' lab to re-compute 
Aa/a for the entire sample. 

There was one complicating aspect of this effect ex- 
cluded from the discussion above, arising from a 7 th spec- 
tral pair. The Sv(X b s )7 showed a more significant non- 
zero slope than the other 6, suggesting a small but signif- 
icant calibration problem with that particular spectrum. 
We therefore applied a slightly more complicated trans- 
formation to the data to allow for this, using a Monte 
Carlo simulation to estimate the potential impact on our 
full combined Keck and VLT sample of both the previ- 
ous effect measured in 6 quasars plus the effect derived 
from the 7 th quasar simultaneously, applied in appropri- 
ate proportions. The full details of this analysis will be 
discussed separately in [7] . 

A systematic of the same magnitude as that from the 
7 th pair cannot be present in any large fraction of our 
data, otherwise it would generate large numbers of no- 
ticeable outliers. If we apply f(Sv) from the 6 quasar 
pairs, the significance of the dipole+monopole model 
Aa/a = AcosG + m, is reduced to 3.1a. Blindly in- 
cluding the effect of the 7 th pair under a Monte Carlo 
method reduces the significance to a most likely value of 



2.2(7. However, in this circumstance we introduce sig- 
nificant extra scatter into the data above that already 
observed, implying that it over-estimates any systematic 
effect of this type. Additionally, the trend of Sv(X b s )i 
against wavelength is different in magnitude and sign for 
each quasar pair, implying that these effects are likely to 
average out for an ensemble of observations. Thus, appli- 
cation of the effect as described above should be regarded 
as extreme in terms of impact on estimating Aa/a. 

Conclusions. — Quasar spectra obtained using 2 sep- 
arate observatories show a spatial variation in the rela- 
tive spacings of absorption lines which could be due to 
an as yet undetected systematic effect, or a dipole vari- 
ation of a. A fit to the dipole gives a significance of 
> 4.2cr, assuming the error bars described above. As- 
suming a dipole interpretation, the two datasets exhibit 
internal consistency and the directions of the indepen- 
dently derived spatial dipoles agree. The magnitudes of 
the apparent Aa/a variation in both datasets also agree. 
A subset of the quasar spectra observed at both obser- 
vatories permits a direct test for systematics. So far, 
none are found which are likely to emulate the appar- 
ent cosmological dipole in a we detect. Consistency with 
other astronomical data is discussed in [T3] . Consistency 
with laboratory, meteorite, and Oklo natural reactor is 
discussed in |15j . Short- wavelength oscillatory variations 
in the wavelength scale such as those reported by [16] , 
|17) . do not significantly impact on our results. To ex- 
plain our results in terms of systematics would require 
at least 2 different, finely tuned, effects. Future similar 
measurements targeting the apparent pole and anti-pole 
directions will maximise detection sensitivity, and fur- 
ther observations duplicated on 2 independent telescopes 
will better constrain systematics. Most importantly, an 
independent technique is required to check these results. 
Qualitatively, our results could violate the equivalence 
principle and infer a very large or infinite universe, within 
which our 'local' Hubble volume represents a tiny frac- 
tion, with correspondingly small variations in the physi- 
cal constants. 
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FIG. 5. Supplementary figure. Same all-sky illustration as in Fig. 1 showing the combined Keck and VLT Aa/a measurements. 
Squares are VLT points. Circles are Keck points. Triangles are quasars observed at both Keck and VLT. Symbol size indicates 
deviation of Aa/a from zero, i.e. Aa/a — ylcos0. The blue dashed line shows the equatorial region. The grey shaded area 
shows the Galactic plane with the Galactic centre indicated as a bulge. More and larger blue squares are seen closer to the 
a-pole (red filled area) and more and larger red circles are seen closer to the a-antipole. 



